Squeezing of light via reflection from a silicon micromechanical resonator 
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We present the measurement of squeezed hght generation using an engineered optomechanical 
system fabricated from a silicon microchip and composed of a micromechanical resonator coupled to a 
nanophotonic cavity. Laser light is used to measure the fluctuations in the position of the mechanical 
resonator at a measurement rate comparable to the free dynamics of the mechanical resonator, and 
greater than its thermal decoherence rate. By approaching the strong continuous measurement 
regime we observe, through homodyne detection, non-trivial modifications of the reflected light's 
vacuum fluctuation spectrum. In spite of the mechanical resonator's highly excited thermal state 
(10,000 phonons), we observe squeezing at the level of 4.5 ± 0.5% below that of shot-noise over 
a few MHz bandwidth around the mechanical resonance frequency of 28 MHz. This squeezing is 
interpreted as an unambiguous quantum signature of radiation pressure shot-noise. 



Monitoring a mechanical object's motion, even with the 
gentle touch of light, fundamentally alters its dynamics. 
The experimental manifestation of this basic principle of 
quantum mechanics, its link to the quantum nature of 
light, and the extension of quantum measurement to the 
macroscopic realm have all received extensive attention 
over the last half century [H |2]. The use of squeezed 
light, with quantum fluctuations below that of the vac- 
uum field, was proposed nearly three decades ago [3 as a 
means for beating the standard quantum limits in preci- 
sion displacement and force measurements. Conversely, it 
has also been proposed that a strong continuous measure- 
ment of a mirror's position with light, may itself give rise 
to squeezed light [4j|5]. In this Letter, we present such a 
continuous position measurement using an engineered op- 
tomechanical system fabricated from a silicon microchip 
and composed of a micromechanical resonator coupled to 
a nanophotonic cavity. Laser light is used to measure the 
fluctuations in the position of the mechanical resonator 
at a measurement rate comparable to the free dynamics 
of the mechanical resonator, and greater than its thermal 
decoherence rate. By approaching the strong continuous 
measurement regime we observe, through homodyne de- 
tection, non-trivial modifications of the reflected light's 
vacuum fluctuation spectrum. In spite of the mechanical 
resonator's highly excited thermal state (10"^ phonons), 
we observe squeezing at the level of 4.5 ±0.5% below that 
of shot-noise over a few MHz bandwidth around the me- 
chanical resonance frequency of cj^/27r ~ 28 MHz. This 
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squeezing is interpreted as an unambiguous quantum sig- 
nature of radiation pressure shot-noise [6 . 

The generation of states of light with fluctuations be- 
low that of vacuum has been of great theoretical inter- 
est since the 1970s [3, 7-9 . Early experimental work 
demonstrated squeezing of a few percent in a large vari- 
ety of different nonlinear systems, such as neutral atoms 
in a cavity [10 , optical fibers [11], and crystals with bulk 
optical nonlinearities ^2 ^3] . This initial research was 
mainly pursued as a strategy to mitigate the effects of 
shot-noise given the possibility of improved optical com- 
munication [7 and better sensitivity in gravitational wave 
detectors [3^^. In recent years, in addition to being in- 
stalled in gravitational wave detectors [14J, squeezed light 
has enhanced metrology in more applied settings [15 . 

The vacuum fluctuations arising from the quantum na- 
ture of light determine our ability to optically resolve 
mechanical motion and set limits on the perturbation 
caused by the act of measurement [16 . A well-suited 
system to experimentally study quantum measurement 
is that of cavity-opt omechanics, where an optical cav- 
ity's resonance frequency can be designed to be sensi- 
tive to the position of a mechanical system. By moni- 
toring the phase and intensity of the reflected light from 
such a cavity, a continuous measurement of mechanical 
displacement can be made. Systems operating on this 
simple premise have been realized in a variety of experi- 
mental settings, such as in large-scale laser gravitational 
wave interferometers [17], microwave circuits with elec- 
tromechanical elements [I8l [19] , mechanical elements in- 
tegrated with or comprising Fabry-Perot cavities j2QH22] , 
and on-chip nanophotonic cavities sensitive to mechani- 
cal deformations [23^ 24 . Conceptually, the same basic 
optomechanical interaction appears when the collective 
motion of an ultracold gas of atoms shifts the resonance 
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frequency of an optical cavity, and recent gas-phase ex- 
periments have shown a variety of optomechanical effects, 
including squeezed light generation [25] . 

The simplest cavity-optomechanical system consists of 
two modes, an optical and a mechanical resonance, and 
is parameterized by their respective resonance frequencies 
(cjo, oJm) and quality factors (Qo, Qm), and the coupling 
rate of intracavity light to the resonant mechanical mo- 
tion (^o)- The Hamiltonian describing the interaction be- 
tween light and mechanics is i^int = hg^a) ax / x^^i^ with 
the amplitude of the mechanical resonator motion being 
given by X = Xzpf • {o^ + where Xzpf = ^Jh|2mQ^iJJ^ is 
the zero-point fluctuation and meff the effective motional 
mass of the resonator. Here a {a^^ and h {p^^ are the 
annihilation (creation) operators of optical and mechan- 
ical excitations, respectively. The optical cavity decay 
rate, k = oJq/Qo^ is the loss rate of photons from the 
cavity and the rate at which optical vacuum fluctuations, 
or shot-noise, is coupled into the optical resonance [26] , 
Similarly, the mechanical damping rate 71 = uJm/Qm is 
the rate at which thermal bath fluctuations couple into 
the mechanical system. In all experimental realizations of 
optomechanics to date, including that presented here, the 
optomechanical coupling rate go has been much smaller 
than the cavity decay rate k,. As such, without a strong 
coherent drive, the interaction of the vacuum fluctuations 
with the mechanics is negligible. 

Under the effect of a coherent laser drive, the cavity is 
populated with a mean intracavity photon number (nc), 
and one considers the optic al flu ctuations about the clas- 
sical steady-state, d \/(^^ + ^- This modifies the 
optomechanical interaction resulting in a linear coupling 
between the fluctuations of the intracavity optical field 
{Xq = d -\- d^) and the position fluctuations of the me- 
chanical system x: i^int = fiGX^x /x^^^. The paramet- 
ric linear coupling occurs at an effective interaction rate 
of G = vvO^o- Through this interaction, the inten- 
sity fluctuations of the vacuum field X^^Q(t) entering the 
cavity impart a force on the mechanical system. 
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which is the radiation pressure shot-noise (RPSN) force. 
The mechanical motion in turn is recorded in the phase 
of the light leaving the cavity. 
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where 9 is the quadrature angle of the detected optical 
field, with = (7r/2) referring to the intensity (phase) 
quadrature. The latter can be interpreted as the infor- 
mational aspect of measurement, with the optical cavity 
playing the role of the position detector (measuring ob- 
servable f at a rate Fmeas = 4G^//^), while the former is 
the quantum measurement back-action imposing a fluc- 
tuating noise force onto the mechanical system [2j . In ad- 
dition to the back-action noise, thermal fluctuations from 



the bath also drive the mechanical motion, with their 
magnitudes becoming comparable as Fmeas approaches 
the thermalization rate FthermaK'^) = 7i^(^)- The bath 
occupation number n{uj) is equivalent to its high temper- 
ature limit k^T})/fux) in these measurements. 

The auto-correlation of the optical output field has con- 
tributions due to the auto-correlation of the optical in- 
put field and the mechanical fluctuations, as well as their 
cross-correlations. Due to the back-action force on the 
mechanical system, the motion has non-zero correlations 
with the RPSN {{x{t)F^x{t')))^ which can be negative 
for some values of 6. Under certain conditions this can 
lead to the squeezing of the optical output field. The 
spectral density of a balanced homodyne measurement of 
the output field, normalized to shot-noise, is given in the 
quasi-static limit for frequencies uo <C uo^i (see the Ap- 
pendix): 



5?f(c^) = l+4(r^eas/a;m)sin(2^) 
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Measurement of the spectrum below shot-noise (squeez- 
ing) is therefore attributable to cross-correlations of mea- 
surement back-action noise due the optical vacuum fluc- 
tuations and position fluctuations of the mechanical res- 
onator. 

The primary technical hurdle to observing such squeez- 
ing, as in many quantum measurements, is the strong cou- 
pling of a preferred detection channel (the optical probe) 
simultaneous with the minimization of unwanted environ- 
mental perturbations of the mechanical system. As indi- 
cated in Eq. ( 11 ), fluctuations from the thermal bath limit 



squeezing of the optical probe field to a regime in which 
n{uo) < Qm- This requirement is equivalent to having a 
Q-frequency product, Qm^m > ^B^b/^- Equation (11) 



also indicates that squeezing becomes appreciable only as 
the measurement rate of the mechanics by the light field 
approaches the mechanical frequency. In some respects, 
both of these challenges have been overcome in recent op- 
tomechanics work [TH [24l : whereby thermal bath cou- 
pling has been made smaller than the optically-induced 
cooling of the mechanical resonator. However, significant 
squeezing over an appreciable spectral bandwidth requires 
not only a large cooperativity between light and mechan- 
ics, as represented by C = Fmeas/Ti and realized in cool- 
ing experiments, but the more stringent requirement that 
the effective measurement back-action force be compa- 
rable to all forces acting on the mechanics, including the 
elastic restoring force of the mechanical structure. A more 
comprehensive model (see Appendix), including the reso- 
nant response of the mechanical system to RPSN, shows 
that the strict requirements of the quasi-static model are 
somewhat relaxed, and that the squeezing scales approx- 
imately as Fmeas/^^ where S is the effective bandwidth of 
squeezing around the mechanical resonance frequency. 

In order to meet the requirements of strong measure- 
ment and efficient detection, we designed a zipper-style 
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[• I( ; . 1: Optomechanical device, a, Scanning electron mi- 
croscope image of a waveguide- coupled zipper optomechanical 
cavity. The waveguide width is adiabatically tapered along its 
length and terminated with a photonic crystal mirror next to 
the cavity. The tapering of the waveguide allows for efficient 
input/output coupling while the photonic crystal termination 
makes the coupling to the cavity single-sided. Two zipper cav- 
ities are coupled above and below the waveguide, each with a 
slightly different optical resonance frequency allowing them to 
be separately addressed, b, (left) Close-up of the coupling re- 
gion between one of the cavities and the waveguide, (right) 
Finite element method (FEM) simulation of the cavity field 
leaking into the waveguide (log scale). Note that the field 
does not leak into the mirror region of the waveguide, c, (top) 
FEM simulation showing the in-plane electrical field of the 
fundamental optical cavity mode, (bottom) FEM simulation 
of the displacement of the fundamental in-plane differential 
mode of the structure with frequency cJm/STT = 28 MHz. The 
mechanical motion, modifying the gap between the beams, 
shifts the optical cavity frequency leading to optomechanical 
coupling. 



optomechanical cavity [23l with a novel integrated waveg- 
uide coupler fabricated from the 220 nm thick silicon de- 
vice layer of a silicon-on-insulator microchip (see Fig.[l^). 
The in-plane differential motion of the two beams at a 
fundamental frequency of uj^/2t: = 28 MHz strongly 
modulates the co-localized fundamental optical resonance 
of the cavity with a theoretical vacuum coupling rate of 
^o/27r = 1 MHz. As shown in Fig. we use a silicon 
waveguide with a high reflectivity photonic crystal end- 
mirror to efficiently excite and collect light from the op- 
tical cavity. Coupling light from the silicon waveguide to 
a single-mode optical fiber is performed using an optical 
fiber taper and a combination of adiabatic mode coupling 
and transformation (see Fig. ^) . 

The experimental setup used to characterize the zipper 
cavity system and measure the optomechanical squeezing 
of light is shown in Fig. |2^. The silicon sample is placed 
in a continuous flow ^He cryostat with a cold finger tem- 
perature of 10 K. A signal laser beam is used to probe 
the optomechanical system and measure the mechanical 



motion of the zipper cavity. A wavelength scan of the re- 
flected signal from the cavity is plotted in Fig.[2]3, showing 
an optical resonance with a linewidth /^/27r = 3.42 GHz 
at a wavelength of Ac = 1540 nm. Inefficiencies in the 
collection and detection of light correspond to additional 
uncorrelated shot-noise in the signal and can reduce the 
squeezing to undetectable levels. For the device studied 
here, the cavity coupling efficiency, corresponding to the 
percentage of photons sent into the cavity which are re- 
fiected, is determined to be r]k = 0.54. The fiber-to-chip 
coupling efficiency is measured at t^cp =0.90. A homo- 
dyne detection scheme [2^ allows for high efficiency de- 
tection of arbitrary quadratures of the optical signal field. 
Characterization and optimization of the efficiency of the 
entire optical signal path and homodyne detection sys- 
tem (see Appendix for details) results in an overall setup 
efficiency of r^setup = 0.48, corresponding to a total signal 
detection efficiency of r^tot = VsetupV^ = 0.26. 

Figure |2]3 shows the noise spectrum of the thermal mo- 
tion of the mechanical resonator obtained by positioning 
the laser frequency near the cavity resonance and tuning 
the relative local oscillator (LO) phase of the homodyne 
detector, ^lock, to measure the quadrature of the refiected 
signal in which mechanical motion is imprinted (roughly 
the phase quadrature for near-resonance probing). The 
mechanical spectrum is seen to contain the in-plane dif- 
ferential mode of interest at a;^/27r = 28 MHz, as well 
as several other more weakly coupled mechanical reso- 
nances of the nanobeams and coupling waveguide (the 
in-plane differential mode peak appears reduced relative 
to the other modes in this plot due to the limited resolu- 
tion bandwidth of the measurement). A high-resolution, 
narrowband spectrum of the in-plane differential mode is 
displayed as an inset to Fig. [2]3, and shows a linewidth of 
7i/27r = 172 Hz, corresponding to a mechanical Q-factor 
of Qm = 1.66 X 10^. The vacuum coupling rate of the 
in-plane differential mode, measured from the detuning 
dependence of the optical spring shift and damping (see 
Appendix), is determined to be ^o/27r = 750 kHz, in good 
correspondence with theory. From the calibration of the 
noise power under the Lorentzian in Fig. |2]3, the in-plane 
differential mode is found to thermalize (at low optical 
probe power) to a temperature of Tb ~ 16 K, correspond- 
ing to a phonon occupancy of (n) ~ 1.2 x 10^. This yields 
a ratio, Qm^^m/^B^b ^ 13, well within the regime where 
coherent motion and squeezing are possible. 

In order to systematically and accurately study the 
noise properties of the reflected optical signal from the 
cavity we make a series of measurements to character- 
ize our laser and detection setup. Figure [2]i shows the 
measured noise power spectral density (PSD) of the bal- 
anced homodyne detector (dark current subtracted) for 
uj ~ ujrn as a function of LO power (signal blocked), 
indicating a linear dependence on power and negligible 
added noise above shot-noise. In the measured squeezing 
data to follow, a LO power of 3 mW is used. Calibra- 
tion of the laser intensity and frequency noise over the 
frequency range of interest {uj /27r = 1-40 MHz) is mea- 
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FIG. 2: Experimental setup and device characterization, a A tunable external- cavity diode laser, actively locked to 
a wavemeter (A-meter), is used to generate a strong local oscillator (LO) and the measurement signal. A portion of each 
field is split off, detected, and used to individually stabilize the intensity in each arm using intensity modulators (IM). Fiber 
polarization controllers (FPC) adjust the polarization of the LO and signal. A variable optical attenuator (VOA) is used to 
set the signal power and an acousto-optic modulator (AOM) is used to generate a tone of known amplitude for calibration (see 
Appendix). The signal beam is fed into a fiber taper that is mounted inside a ^He continuous flow cryostat and is coupled to 
the zipper cavity. The reflected signal beam from the cavity is separated using a circulator, and is switched between one of 
three detection paths: a power meter (PM) for power calibration, a photodetector (PDl) for optical spectroscopy of the cavity, 
and an erbium-doped fiber amplifier (EDFA) for spectroscopy of the coupled cavity-optomechanical system on a real-time 
spectrum analyzer (RSA) or network analyzer (ENA). For the measurement of squeezing, the reflected cavity signal is sent to 
a fourth path containing a variable coupler (VC), where it is recombined with the strong local oscillator and detected on a 
balanced homodyne detector (BHD). The relative phase between the LO and the signal is set using a fiber stretcher (FS). b 
(top) Reflected signal from the optical cavity at low optical power (ric ~ 10, linewidth /^/27r = 3.42 GHz), (bottom) High- 
power (ric = 790) reflected signal, showing the cavity-laser detuning (dashed line) locked to during squeezing measurements, c 
Homodyne noise PSD of the reflected signal showing the transduced thermal Brownian motion of the zipper cavity at Tb = 16 K 
(green curve; ric = 80). The red curve is the shot noise level and the black curve is the detector's dark noise. The inset shows 
a zoom-in of the fundamental in-plane differential mechanical mode of the zipper cavity (linewidth 7i/27r = 172 Hz), d Mean 
value of the PSD of the balanced homodyne detector as a function of the LO power (signal blocked). The filled data point 
indicates the LO power used in the squeezing measurements. The red and dashed black curves correspond to a linear fit to 
the data and the detector dark current level, respectively, e Noise PSD as a function of ^lock with the signal detuned far 
off-resonance at A/k ^ 30 (shades green to red), referenced to the noise level with the signal blocked (blue). 



sured by direct photodetection of the laser, pre- and post- 
transmission through a fiber Mach-Zehnder interferome- 
ter with a known frequency response. The laser intensity 
noise is measured to be shot-noise dominated over this fre- 
quency range, while the laser frequency noise is measured 
to be roughly flat at a level of S^ju ~ 5 x 10^ rad^-Hz. 
Laser phase noise on the signal beam can be converted 
to intensity noise by reflection from the dispersive cavity 
or due to frequency dependent components in the optical 
train, and can add to the detected noise floor. Due to 
the broad linewidth of the zipper cavity resonance, the 
effects of the laser phase noise are found to be negligi- 
ble. The full suite of noise and calibration measurements 



performed are described in detail in the Appendix. 

Measurement of the noise in the reflected optical signal 
from the cavity as a function of quadrature angle, fre- 
quency, and signal power is presented in Figs. [3] and [4] 
These measurements are performed for laser light on res- 
onance with the optical cavity and for input signal powers 
varying from 252 nW to 3.99 /iW in steps of 2 dB, with the 
maximum signal power corresponding to an average intra- 
cavity photon number of (ric) = 3, 153. Positioning of the 
laser at the appropriate cavity detuning for each signal 
power is performed by scanning the wavelength across 
the cavity resonance while monitoring the reflection, and 
then stepping the laser frequency towards the cavity from 
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Q|oc/^ ^loc/^ Frequency (MHz) 

FIG. 3: Optomechanical squeezing of light, a, Theoretical model. Density plot of the predicted reflected signal noise 
PSD, as measured on a balanced homodyne detector and normalized to shot-noise, for a simplified model of the optomechanical 
system (see Appendix). Areas below shot-noise are shown in blue shades on a linear scale. Areas with noise above shot-noise 
are shown in orange shades on a log-scale. The solid white line is a contour delineating noise above and below shot-noise, b, 
Experimental data. Density plot of the measured reflected signal noise PSD for ric = 790 normalized to the measured shot-noise 
level, c. Slice of the measured density plot in b taken at ^lock/vr = 0.23. d. Slice of the measured density plot in b taken at 
Oiock/TT = 0.16. In c and d, the black curve corresponds to the measured data slice extracted from b. The dark blue traces are 
several measurements of the shot- noise level (average shown in light blue). Also indicated is a model of the squeezing in the 
absence of thermal noise (red dashed curve), the same model with thermal noise included (solid red curve), and a full noise 
model including additional phenomenological noise sources (solid green curve). 



the red side until the reflection matches the level that cor- 
responds to a detuning of Aiock/^ ^ 0.04 (see Fig. 
The laser is locked to this frequency using a wavemeter 
with a frequency resolution of ±0.0015/^. Drift of the op- 
tical cavity resonance over the timescale of a single noise 
spectrum measurement (minutes) is found to be negligi- 
ble. An estimate of the variance of Aiock is determined 
from the dependence of the transduction of the mechani- 
cal motion on the quadrature phase, indicating that from 
one lock to another Aiock//^ = 0.044 ± 0.006. 

In Fig. |3] we plot the theoretically predicted and mea- 
sured noise PSD versus quadrature angle for a signal 
power corresponding to (ric) = 790 photons. Each 
quadrature spectrum is the average of 150 traces taken 
over 20 seconds, and after every other spectrum, the sig- 
nal arm is blocked and the shot-noise PSD is measured. 
The shot-noise level, which represents the noise of the 
electromagnetic vacuum on the signal arm, is used to nor- 
malize the spectra. We find at certain quadrature angles, 
and for frequencies a few MHz around the mechanical res- 
onance frequency, that the light reflected from the zipper 
cavity shows a noise PSD below that of vacuum. The den- 
sity plot of the theoretically predicted noise PSD (Fig.|3^) 
shows the expected wideband squeezing due to the strong 
optomechanical coupling in these devices, as well as a 
change in the phase angle where squeezing is observed at 
below and above the mechanical frequency. This change is 
due to the change in sign of the mechanical susceptibility 



and the corresponding change in phase of the mechanical 
response to RPSN. The measured noise PSD density plot 
(Fig. [sj)) shows the presence of several other mechanical 
noise peaks and a reduced squeezing bandwidth, yet the 
overall phase- and frequency-dependent characteristics of 
the squeezing around the strongly-coupled in-plane me- 
chanical mode are clearly present. In particular. Figs. [Sj^ 
and d show two slices of the noise PSD density plot which 
show the region of squeezing change from being below to 
above the mechanical resonance frequency. 

In Fig. [4^ we show the measured noise PSD (grey cir- 
cles) as a function of quadrature angle for a frequency 
slice at cj/27r = 27.9 MHz of the data shown in Fig. [3J3. 
The measured squeezing (anti-squeezing) is seen to be 
smaller (larger) than expected from a model of the op- 
tomechanical cavity without thermal noise. We also plot 
in Fig. the maximum measured and modeled squeez- 
ing as a function of signal power. The simple theory pre- 
dicts a squeezing level (blue curve) which monotonically 
increases with signal power, whereas the measured max- 
imum squeezing saturates at a level of 4.5 ± 0.5% below 
the shot-noise at an intracavity power (ric) = 1,984 pho- 
tons. In order to understand the noise processes that limit 
the bandwidth and magnitude of the measured squeez- 
ing, we plot in Fig. ^ the noise PSD (shot-noise sub- 
tracted) for phase quadratures that maximize (left plot) 
and minimize (right plot) the transduction of the me- 
chanical mode peak at a;^/27r = 28 MHz. Along with 
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FIG. 4: Spectral and power dependence of noise, a, 

Measured (filled circles) balanced homodyne noise power of 
the reflected signal at uj /27v = 27.9 MHz versus quadrature 
angle (Aiock/^ = 0.044 and ric = 790). The red curve cor- 
responds to the full noise model. The solid blue curve is for 
a model including the response of the mechanical mode in 
the absence of thermal noise, i.e., driven by RPSN only (the 
dashed blue curve shows the thermal noise component), b. 
Measured (filled circles) minimum noise PSD normalized to 
shot-noise versus ric. The left plot is the maximum squeezing 
for (jj < ujm and the right for uj > ujm- Also shown is the 
single-mode noise model (blue curve) and the full noise model 
(red curve), c. Balanced homodyne noise PSD of the reflected 
cavity signal for Aiock/^ = 0.052 and Uc — 3, 153. Left (right) 
plot shows phase quadrature corresponding to maximum (min- 
imum) transduction of mechanical motion. The black curve is 
the measured data with the shot-noise level subtracted. Also 
shown are modeled laser phase noise (green curve), the single- 
mode noise model (blue curve), and the full noise model (red 
curve) . 



the measured data (black curve), we also plot the esti- 
mated noise due to phase noise of the signal laser (green 
curve) and that for a single mechanical mode (blue curve) 
assuming a thermal bath temperature of T5 = 16 K and 
a frequency-independent damping rate. The single- mode 
noise model greatly underestimates the background noise 
level, especially in the quadrature minimizing transduced 
motion in which we measure a cj"^/^ (as opposed to o;"^) 



frequency dependence to the low frequency noise. As de- 
scribed in the Appendix, the additional measured noise is 
thought to arise from a combination of the thermal noise 
tails of higher frequency mechanical modes and fluctu- 
ations in the optical cavity damping rate, along with a 
small amount of heating due to optical absorption in the 
silicon cavity. The red curves in each of the plots in Fig. [4] 
show the full noise model incorporating these added phe- 
nomenological terms. 

Somewhat surprisingly, these measurements show that 
by reflecting light off a thin-film mechanical resonator un- 
dergoing large amplitude thermal motion, light that is in 
certain respects quieter than vacuum can be obtained. 
This is found to result from radiation-pressure fluctua- 
tions being strongly imprinted on and modified by the 
motion of the mechanical resonator, thus demonstrating 
a fundamentally quantum aspect of displacement mea- 
surement. The devices in this work utilize lithographic 
patterning at the nanoscale to transform silicon into a 
material with an effective quantum optical nonlinearity at 
an engineerable optical wavelength. The modest level of 
squeezing realized in this work is predominantly limited 
by thermal noise, but also by the efficiency with which 
the reflected light can be collected by external optics. 
The effects of thermal noise may be substantially reduced 
by working with materials of higher intrinsic Q-frequency 
product, such as diamond and silicon carbide [29 . Given 
the microchip form of the devices studied here, and the 
potential for device integration, it is interesting to con- 
sider whether new squeezed light applications might arise. 
For example, squeezed light generated by one device could 
be directly sent into another device for use as an optical 
probe. Such an on-chip squeezer and detector could be 
used as a quantum-enhanced micro-mechanical displace- 
ment and force sensor [30 . More generally, we expect 
future experiments with feedback and strong measure- 
ment of the dynamics of a mechanical system to be within 
reach. In addition, using quantum light as an input we 
expect to be able to generate entangled states of mechan- 
ics and light with similar devices. 
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I. THEORY 

Optomechanical systems can be described theoretically 
with the Hamiltonian (see main text) 

H = hWod)a + humob^b + hgoa^a{b^ + 6), (Al) 

where d and b are the annihilation operators for pho- 
tons and phonons in the system, respectively Generally 
the system is driven by intense laser radiation at a fre- 
quency cjl , making it convenient to work in an interaction 
frame where Uq is replaced by A in the above Hamilto- 
nian with A = — ^L- To quantum mechanically de- 
scribe the dissipation and noise from the environment, we 



(iA+^)a-i^oa(6t + 6) 



igod^d - v^6in(t). 



which account for coupling to the bath with dissipation 
rates /^i, and 71 for the intrinsic cavity energy decay 
rate, optical losses to the waveguide coupler, and total 
mechanical losses, respectively The total optical losses 
are = Kjq ~\~ Kj\ . These loss rates are necessarily accom- 
panied by random fluctuating inputs ain(t), ain,i(t), and 
bin{t)^ for optical vacuum noise coming from the coupler, 
optical vacuum noise coming from other optical loss chan- 
nels, and mechanical noise (including thermal). 

The study of squeezing is a study of noise propagation 
in the system of interest and as such, a detailed under- 
standing of the noise properties is required. The equa- 
tions above are derived by making certain assumptions 
about the noise, and are generally true for the case of 
an optical cavity, where thermal noise is not present, and 
where we are interested only in a bandwidth of roughly 
10^ smaller than the optical frequency (0 - 40 MHz band- 
width of a 200 THz resonator). For the mechanical sys- 
tem, where we operate at very large thermal bath oc- 
cupancies (> 10^) and are interested in the broadband 
properties of noise sources (0 - 40 MHz for a 30 MHz 
resonator), a more detailed understanding of the bath is 
required, and will be presented in the section on thermal 
noise. 

At this point, we linearize the equations assuming a 
strong coherent drive field ao, and displace the annihi- 
lation operator for the photons by making the transfor- 
mation d ao -\- d. This approximation, which neglects 
terms of order is valid for systems such as ours where 
go <C i.e. the vacuum weak coupling regime. We are 
then left with a parametrically enhanced coupling rate 
G = ^o|<^o|- Using the relations given in the mathemati- 
cal definitions section (VI) of this document, we write the 



solution to the QLEs in the Fourier domain as 



d{uj) 
b{u) 



/^ain(cj) - y^ain,i(cj) - iGibjuj) + b\uj)) 
7[6in(a;) iG{d{uj) + at (a;)) ^^^^ 



z(cjmO -^) +7i/2 ^(^mO " ^) +7i/2* 



Finally we note that by manipulation of these equa- 
tions, the mechanical motion can be expressed as a (renor- 
malized) response to the environmental noise and the op- 
tical vacuum fluctuations incident on the optical cavity 
through the optomechanical coupling 
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i{uJm -^) +7/2 

iG y^ain(cj) + y^ain,i(cj) 

i{A - cj) + i{ujm -oj) ^ 7/2 

iG y/i^al^juj) + V^aL,i(^) 

-z(A + cj) + /^/2 i(cJm-^)+7/2 



(A3) 



The renormalized mechanical frequency and loss rate 
are uj^i = ^mo + ^^m, and 7 = 7i + 70M, respectively, 
with 



6ujm = IGrim 



70M = 2|GpRe 



1 



1 



i{A - UJm) + ^/2 -i(A + UJm) + /^/2 

(A4) 

1 1 



i(A - CJm) + ^/2 -i(A + CJm) + 

(A5) 



It is convenient to define here what we mean by a 
quadrature, as it is the observable of the light field that 
our measurement device (the balanced homodyne detec- 
tor (BHD) setup) is sensitive to: 



(A6) 



We are interested in the properties of for various 

quadrature angles ^, given the influence of the mechanical 
system. 

The measurement of the field provides us with a record 
i(t) = XQ^^^\t) for a certain 6. We use a spectrum an- 
alyzer to perform Fourier analysis on this signal and ob- 
tain a symmetrized classical power spectral density (PSD) 
Sn{u), as defined in the mathematical appendix (section 

For a vacuum field such as the input field, the measured 
quadrature X^^^^^ (t) will have a power spectral density 



1. 



(A7) 



This is the shot-noise level which is due to the quan- 
tum fiuctuations of the electromagnetic field. Mathemat- 
ically, it arises from the correlator (avac(^)^tac('^O) — 
S{uj + a;'), with all other correlators {al^^{uj)a^s,c{^^)) ^ 
{al^^{uj)al^^{uj')), (avac(^)^vac(^0)5 arising in the expres- 
sion {r {uj)I{uj')) equal to zero. 



A. Approximate quasi-static theory 

In this section we present a simplified derivation of how 
squeezing is obtained in the studied optomechanical sys- 
tem to elucidate the important system parameters and 
their role in squeezing. We make a few approximations 
to simplify the derivation: 



1. A = 0: The laser is tuned exactly to the optical 
cavity frequency. 

Perfect coupling. 

3. hz:^ uo^: Bad cavity limit. 

4. cj <C cjm- We are only interested in the quasi-static 
response, so the resonant response of the mechanical 
resonator does not play a role. 

Under these assumptions, equations ([2| and (|3| can be 
written as (using the relation for the optical output field 
ttout(^) = ^in(^) + V^a(a;)): 

iujj){uj) = -V7i6in(a;) + ^(ain(a;) + aL(^))' 

aout(w) = -ai„(w) - '^{b{uj) + b\uj)). (A8) 

The first equation shows the response of the mechani- 
cal resonator subsystem to the thermal bath fiuctuations 
(^in(^)) and the optical vacuum noise from the measure- 
ment back-action. We define Fmeas = / and inter- 
pret it as the measurement rate [2], such that the factor 
appearing in front of the optical vacuum noise operators 
is vTmeas- This rate also appears in the second equation 
for the output field, in front of the normalized position 
operator x/x^^f = b{uj) + b^{uj), which is the observable 
that is being measured. 

Note, from the expression for aout(^) it follows, that 
since the position is a real observable with an imaginary 
prefactor, the effects we consider depend strongly on the 
quadrature being probed, i.e. the real part of the expres- 



sion, X, 



(out) 
=0 ' 



will not be affected by the optomechanical 



coupling. 

At this point we can easily calculate the properties of 
the detected spectrum Sjj^{uj)^ by writing aout in terms 
of ttin and 6in for which the correlators are known: 

^out(^) = -^in(^) - 2i^^^^^(ain(a;) + 



+ 



'-{bM-bliu)). (A9) 



Ignoring thermal noise for the moment (71 = 0), and drop- 
ping terms of order (Fmeas/^m)^ (assuming Fmeas < ^m) 

we arrive at: 

/oo 
-00 

= i + 4(r 

meas 

/cjm)sin(2(9). (AlO) 

Note that for certain values of 6>, the detected spectral 
density can be smaller than what one would expect for a 
vacuum field. For = — 7r/4, we achieve the maximum 
squeezing with a noise fioor of 1 — 4(Fmeas/^m) which 
strongly dependends on the ratio Fmeas/^m- 

To understand the effect of thermal noise, we as- 
sume the form of the correlator to be (6in(co')6-^(co'')) = 
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Measurement Phase (ti radians) Frequency (MHz) Frequency (MHz) 



Figure Al: Squeezing theory, a, Density plot of the predicted squeezing SjY^{uj) vs. phase angle and frequency normalized 
to the shot-noise. The mechanical mode can clearly be seen at cJm/STr = 30 MHz. The solid white lines outline the region 
where the power spectral density falls below 1 (the shot-noise level) indicating the presence of squeezing for that phase and 
frequency. The dashed white lines at ^ = — 7r/4 and = +7r/4 correspond to regions where squeezing can be obtained below 
and above the mechanical frequency respectively and the components of the noise model for these phases is shown in detail 
in figures b and c. In these figures the spectra are again normalized to the shot-noise level plotted as a grey line. The simple 
squeezing model without thermal noise (Eq. (^^[lO|) is represented by the dashed green line and the simple model with thermal 
noise (Eq. (A[TTt) is the sohd green line. The solid black line is the full squeezing model SjY^{(jj) corresponding to a with 
the constituent components: the contribution from the optical vacuum fluctuations {Sj^^ai^)'^ Eq. (A |24| )) represented by the 
dashed black line and the thermal noise (^^//^^(cc;); Eq. (i^[25t) represented by the dashed red line. 



{n{uj) + l)^(a; + u'), {bl^{uj)kn{c^')) = n{uj)d{uj + u'), 

{Hni^)Hni^^)) ~ 0' (^in(^ )^in(^ O) = (these cxprcs- 
sions are discussed in section IV G). Then a calculation 



similar to the one leading to equation ( |10| ) gives 

S^fiuj) = l+4(r^eas/c^m)sin(2^) 



+4- 



. n{uj) 



(l + cos(2^)), (All) 



where we have assumed n(co'), the bath occupation at fre- 
quency cj, to be much larger than unity. At = — 7r/4, 
we have 

(C^) = 1 - 4(r^eas/c^m)(l - n{uj)/Qm)^ (Al2) 

In this model, there is no squeezing at u if n{uj) > Qm- 
This means that for squeezing to be present, coherent 
evolution of the mechanical system must be possible, i.e. 
the rate at which phonons enter the mechanical system 
from the bath (71 n) must be smaller than the mechanical 
frequency uo^. In conclusion, the important requirements 
to achieve squeezing are to make Fmeas comparable to 
cjm, and to reduce the thermal occupancy or increase the 
mechanical quality factor to achieve n{uj) < Qm- 



B. The effect of dynamics and correlation between 
RPSN and position 

As a next step, we take into account the dynamics of the 
mechanical resonator while keeping the approximations of 
the bad-cavity limit {hz ujrn) and on-resonant probing 
(A = 0). In addition to further clarifying some of the 
observed features, this treatement, as presented in the 



main text, elucidates the role of correlations between the 
mechanical system's position and the back-action force. 

The response of the mechanical system to a force is 
captured by its susceptibility: 



Xm(^) 



1 



(A13) 



The form of the damping considered here is the strongly 
sub-ohmic structural damping which is observed in our 
measurements [33l|34] (cf. Section IV G). The mechanical 
system responds to random noise forces F^Yit) from the 
thermal bath (which we treated in the last section and 
neglect here), and to the quantum back-action from the 
cavity Fba(^). 

The back-action force for the resonant case can simply 
be found by linearizing the expression for the radiation 
pressure force FRp(t) = —hga^a/xzpf. We find the force 
imparted on the mechanics due to the shot-noise of the 
cavity field to be 



FBA{t) = 



meas V>(in) 



(A14) 



for the case of resonant driving. The fiuctuations im- 
parted on the mechanics are from the intensity quadra- 
ture of the light (^ = 0). Using equation (|8|, we can write 
the output field quadrature as: 



r(in) / 



^zpf 



-f (t) •sin(i9). (A15) 



We note here that the mechanical position fiuctuations 
are primarily imprinted on the phase quadrature of the 
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output light, with 6 = ±7t/2. The intensity quadrature 
is unmodified {Xg^Q\t) = — X^^^Q(t)) since changes in 
the cavity frequency are not transduced as changes in 
intensity when the laser is resonant with the cavity. 
The output of the homodyne detector normalized to the 



shot-noise level is found by taking the auto-correlation of 
eqn. (A 15). The correlations between radiation pressure 
shot-noise and the mechanical motion are important in 
this calculation [35-41 and must be taken into account. 
In the time-domain we find the auto-correlation to be: 



-^zpf 

+2^-^ sin(i9) cos{0){FBA{t)x{f) + x{t)FBA{t')). (A16) 



The cos(^) in the last term comes from the general ex- 
pression for a quadrature XQ^\t) = XQ^Q{t) cos{6) + 
XQ^]^^2i^) ^'^^i^)^ equation (A[14|. The key compo- 
nents of equation (A 16 ) are the shot-noise level, the ther- 
mal noise, and the cross-correlation between the back- 
action noise force and mechanical position fiuctuations. 
It is only the latter which can give rise to squeezing, 
by reducing the fiuctuation level below shot-noise. This 
squeezing can be calculated spectrally: 



5'sq(^) 



/CO 
dr[{FBA{t)x{t-r)) 
-oo 



+ {x{t)FBA{t - T))]e«^- 
= 2;isin(2^)%^^XmM 

^zpf 

= 4mWmSin(20)rmeasXm(w). 



(A17) 



At the DC or quasi-static limit {uj 0) the susceptibility 
Xm can be used and we reobtain the results 



IA| (cf. equation (A[To|). We see that for 



from section 

6> < 0, squeezing is obtained in this limit. At frequencies 
larger than c^m, Xm('^) changes sign, and we expect to see 
squeezing at quadrature angles > 0. Additionally, since 
Xm(^) becomes larger around the mechanical frequency, 
we expect the maximum squeezing to be enhanced. More 
specifically, at a detuning S = ujrn—^ from the mechanical 
resonance, we expect the parameter characterizing the 
squeezing to be proportional to Fmeas/^^- All of these 
features are seen in Fig. [l] 

It is important to note here that in the absence of other 
nonlinearities in the system, any reduction of the noise 
below the vacuum fiuctuations can only be caused by the 
correlations between the RPSN and the position fiuctua- 
tions of the system. This makes the problem of proving 
the correlations between RPSN and mechanical motion 
equivalent to the problem of proving that the refiected 



light from the optomechanical cavity has been squeezed. 

Conceptually this form of probing the RPSN is similar 
to that carried out by Safavi-Naeini et al. [38", "^T and ana- 
lyzed by Khalili et al. [39 . It also shares features with the 
cross-correlation measurements proposed by Heidmann et 
al. [35 , and B0rkje et al. [37 , and recent experiments by 
Purdy et al. [6]. The distinguishing feature of this type 
of measurement is that the quantum correlations between 
the fiuctuations of the position and the electromagnetic 
vacuum manifest themselves as squeezed light. 



C. General derivation of squeezing 

Among the approximations made in section |I A[ the 
quasi- static approximation is the least correct. In fact, 
in our experiments, the most observable squeezing occurs 
with UJ close to c^m and even slightly larger than cjm, so 
OJ <C oJm is not valid. Near the mechanical frequency, 
resonant enhancement of the optical vacuum fiuctuations 
by the mechanical resonator causes squeezing greater than 
that predicted in the quasi-static regime to be possible. 

Here we show the results of a derivation that does not 
rely on most of the assumptions used in the approximate 
model. Of the assumptions in the previous section, the 
only simplification we keep here is to assume perfect cou- 
pling The effect of imperfect coupling can be 
taken into account trivially and is explained after this 
section (see I C 1 ). 

By substitution of equation (|3| into the equation for 
a{uj) ([2|, we arrive at: 



-^B,{uj)bM + B2{uj)bl{uj), (A18) 



with 
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i{A-uj) + /^/2 



1 



_i(A -Lo)-\- tv/2 i{Lo^ -co)-\- 7/2 i(A - Cl;) + a^/2 -^(cJm + c^;) + 7/2 



i(A - cj) + ^/2 



-i(A + dj) + /^/2 i{u;^ - cu) -\- 7/2 -i(A + dj) + k/2 + Cc;) + 7/2 



//^7i 



i(A-cj) + /^/2 



i(cJni -^) +7/2. 

iG 



i(A - (j) + tv/2 [-i{cOm -\-co) -\- 7/2_ 



(A19) 

(A20) 
(A21) 
(A22) 



The expressions give us the output field in terms of the We can calculate Sjj^{uj) from this expression, which we 
input fields, since split into two parts, one only due to the optical vacuum 

fluctuations, and the other containing the contribution 
aout(^) = ain(^) + v^a((j)) from thermal noise: S'jfiuj) = S'^yU^) + 5'?y,U^)- 

^Bi{uj)kn{uj) + B2{uj)bl{uj). (A23) 



S^jj^icj) = \A2{-uj)\' + |1 + AiMp + 2Re{e-2^^(l + A,{uj))A2{-cu)} (A24) 
S-jj^iuj) = |5iHp(n(a;) + l) + |5i(-a;)pnM 

-^\B2{-oj)\\n{uj) + 1) + \B2{uj)\^n{uj) 

^2Re{e-^'^Bi{uj)B2{-uj)}{n{uj) + 1) + 2Re{e-^'^ Bi{-uj)B2{uj)}n{uj) (A25) 



1. The effect of imperfect optical coupling and inefficient 
detection 

At every juncture in an experiment where the opti- 
cal transmission efficiency is less than unity (77 < 1), 
an equivalent optical circuit can be defined involving an 
77 : (1 — 77) beam splitter with the output being r] times the 
input and (1 — 77) times the vacuum. Therefore the effect 
of optical losses and coupling inefficiencies on the detected 
spectra can be calculated by replacing the measured field 
quadrature with: 

= + yr^X^""^^ (A26) 

This source of vacuum noise is completely unrelated to the 
cavity output, and there are no cross-correlation terms, 
so the detected current spectral density will be given by 

Sfficv) = rjSjTi^) + (1 - v)S}r{^), (A27) 

where S}j^{u;) = 1 is the shot-noise. 

Measurement inefficiencies take two forms, one is due to 
ineffeciencies in the detection, while the second is because 



of excess electronic noise or "dark current" present due 
to the circuitry of the detector and amplifier. This excess 
noise can also be thought of as a detection inefficiency 
by considering the amount of optical shot-noise inserted 
into the signal which would produce it. Since the dark- 
current is measured with no optical input, and the real 
shot-noise level increases linearly with the local oscillator 
(LO) power, this inefficiency is power dependent and can 
be minimized for large LO powers. In our case, the dark 
current was found to be 10.4 dB below the detected shot- 
noise. The total detector efficiency was measured to be 
77HD = 66%. 



II. EXPERIMENT 

A. Measurement of losses 

In order to estimate the total squeezing expected in 
our setup we carefully characterize all losses in our sys- 
tem. Some of these losses are static (e.g. circulator losses) 
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Figure A2: Experimental setup. A detailed description of the experimental setup can be found in the main text. 



while others can vary from experiment to experiment (e.g. 
coupling efficiency of the fiber taper to the waveguide). 
In figure [5] typical losses are shown as efficiencies {r]) for 
various parts of the experiment. The efficiency of send- 
ing light from port 1 to 2 of our optical circulator is 
r]i2 = 85%, and r]23 = 88% for port 2 to 3. In addition, 
the efficiency from port 3 of the circulator to the homo- 
dyne detector is tjsyl = 92%. All these losses are fixed and 
do not change over time as the components are optically 
spliced together. Measuring the coupling efficiency of the 
fiber taper to the waveguide is done every time a new data 
set is taken. This is accomplished by switching the light 
that is reflected from the waveguide to a power meter and 
comparing the reflected power to a known input power 
with the laser tuned off-resonance from the optical mode 
(off-resonance the device acts as a near-perfect mirror). 
Typical achieved efficiencies are around t^cp = 90%. The 
efficiency of the homodyne detection strongly depends on 
the alignment of the polarization between the local oscil- 
lator and the signal, as well as by how much the power in 
the LO overcomes the electronic noise floor of the detec- 
tor. To determine this efficiency we use an acousto-optic 
modulator (AOM) inserted in our setup before the circu- 
lator in the signal path. The AOM shifts the frequency 
of the light creating a tone 88 MHz away from the signal 
with a fixed, known amplitude, and identical polariza- 
tion to the signal (we directly measure the power of this 
tone with a power meter). This tone can now be used 
to determine the total homodyne efficiency by measur- 
ing its power on the spectrum analyzer, taking the other 
losses into account. Our typical homodyne efficiency is 
^HD = 66% resulting in a total setup efficiency (detection 
efficiency of optical signal photons in the on-chip waveg- 
uide) of roughly r/setup = ^CP • V23 ■ Vsu ■ ^hd ^ 48%. 



B. Data collection procedure 

Careful calibration of our data is crucial in understand- 
ing all noise sources and potential drifts over time in our 
setup. The losses in our setup are determined before we 
make a new data run as described in the previous section. 
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Figure A3: Phase information. The blue trace shows the 
interference signal of the local oscillator and the signal on the 
homodyne detector when their relative phase is scanned using 
a fiber stretcher in the local oscillator (LO) arm. The voltage 
reading here is proportional to cos(^ — 0) where is the phase 
difference between the LO and input to the cavity, and is 
the phase imparted by reflection off the cavity. This inter- 
ference signal is used to actively stabilize the relative phase 
to different set points (green traces). Occasionally the lock 
fails, as shown by the red traces, and any associated data is 
discarded. The range in which the phase can be stably locked 
is slightly smaller than —7v/2 to 7r/2 due to the turning points 
in the sinusoidal interference curve. 



We then proceed to record an optical trace of the cavity 
resonance by switching the light to a photodetector (PDl 
in figure [2| and scanning the laser wavelength. This trace 
provides the information to lock the laser to a fixed detun- 
ing (typically 0.04 • hi red of the cavity resonance), which 
is accomplished using a simple software lock and feed- 
back from the wavemeter (with a resolution of roughly 
0.003 • hz) and is described in more detail in the subsec- 
tion below. As a next step the optical signal is switched 
to the homodyne detector and the relative phase between 
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the signal and the local oscillator is scanned using the 
fiber stretcher in the LO arm. The resulting interference 
is shown in figure [3] as the blue trace. The interference 
signal is used to lock the relative phase between the signal 
and LO using a Toptica DigiLock 110. The green traces 
show the properly locked signal, while the red traces are 
phase set points where the lock failed requiring the asso- 
ciated data to be discarded. We then record the spectra 
of the homodyne signal and for every trace taken we also 
save a spectrum of the shot-noise by switching the signal 
arm away from the homodyne detector and only measure 
vacuum input to the signal arm of our detector. We re- 
lock the laser with respect to the cavity every other data 
point to counteract drift. This procedure is repeated for 
several different phases and different input powers. We 
typically took data for 60 different phases for every input 
power within a range of a little less than — 7r/2 to 7r/2. 



C. Relation between detuning and quadrature 

The laser frequency is positioned at a detuning of 
roughly 0.04- by starting at a larger detuning on the red 
side of the cavity, and stepping the laser blue in 0.1 pm 
steps (12 MHz) towards the cavity while monitoring the 
average intensity of the reflected light on PDl. Once 
the target intensity is reached, the laser is kept at this 
wavelength during the course of the measurement by the 
wavemeter lock without further feedback from PDl. The 
intensity reading gives us an idea of the value of the de- 
tuning which is determined more accurately by analysis 
of our homodyne spectra. 

The homodyne spectra are taken at different phase lock 
points (see Fig. |3| corresponding to quadrature angles 
^lock between the reflected signal and the LO. These an- 
gles differ from our convention in Section [T] where the 
phase 6 between the input light into the cavity and the 
local oscillator is considered. They are related to onean- 
other by the phase imparted on the input light upon re- 
flection from the cavity. 



(/)(A) = Arg 
and the relation 



1 



iA + k/2 



(A28) 



(A29) 



For a given laser-cavity detuning A, we sweep through 
the different phase lock points (see Fig. 3| 6>iock, and 
take mechanical spectra for each phase. The phase that 
minimizes the mechanical signal 0^^^^ is determined from 
the recorded spectra. This allows us to solve for A us- 
ing the expression 0^^^^ = r(A) - 0(A) where i9*(A) 
is the phase minimizing the mechanical transduction ac- 
cording to the model in the previous section. To first 
order (for A <C a>:) ^* is since no mechanical signal 
is observed in the intensity quadrature of the reflected 
light. This post-processing of the data allows us to deter- 
mine that across the measured powers the detuning was 
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Figure A4: Detuning and phase lock points, a, An optical 
scan taken before the data run starts is shown. The blue ver- 
tical line denotes the target detuning the software lock moves 
the laser to, determined from the measured reflection inten- 
sity. The laser is kept at that detuning via a wavemeter lock, 
as the light is switched away from PDl and to the homodyne 
detector. The measured area under the mechanical mode is 
plotted in b (blue circles) at this detuning. A minimum value 
is reached for a local-oscillator to reflected signal phase of ^fock- 
Depending on the detuning, different mechanical mode ampli- 
tudes can be measured at this phase angle ^lock, according to 
the model. We obtain an accurate estimate of the detuning by 
calculating the detuning at which the mechanical mode am- 
plitude is minimized at the measured ftock as shown in c. The 
expected mode amplitudes for the detunings represented by 
the red and green lines in c are shown by the similarly colored 
curves in b. 



A = 0.044 • ± 0.006 • tv. For a single measured power, 
we expect a more accurate determination, with an uncer- 
tainty on the order of 0.003 • K:. This level of accuracy in 
the detuning also determines the uncertainty in quadra- 
ture angle of 0.04 rad. 



III. SAMPLE FABRICATION AND 
CHARACTERIZATION 

A. Fabrication 

The devices are fabricated from a silicon on insulator 
(SOI) wafer (SOITEC, 220 nm device layer, 3 /im buried 
oxide, device layer resistivity 4 — 20 Q ■ cm) using elec- 
tron beam lithography and subsequent reactive ion etch- 
ing (RIE/ICP) to form the structures. The buried oxide 
is then removed in hydrofluoric acid (49% aqueous HF so- 
lution) and the devices are cleaned in a piranha solution 
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(3:1 H2SO4 and H2O2) and finally hydrogen terminated 
in diluted HF. For a more detailed description see [42 . 

B. Optical Characterization 

The optical characterization of our devices is done 
by sweeping the laser frequency across the optical res- 
onance while detecting the reflected light in a photode- 
tector (PDl in figure [2|. This light is simultaneously sent 
to a wavemeter to record the absolute wavelength and ac- 
curately determine the linewidth and center frequency of 
the resonance. Each chip contained several designs where 
the waveguide loading (coupling) of the optical cavity was 
varied by changing the gap size between the waveguide 
and nanobeam. For our measurements we chose a slightly 
overcoupled (/^e/^i ^ 1-22 > 1) device with good optical 
quality (57,000 loaded Q) gj. 

C. Mechanical Characterization 



The intrinsic mechanical damping rate 71 and the op- 
tomechanical coupling rate go are measured by detecting 
the mechanical response to the signal laser, through the 
reflected signal field, on the spectrum analyzer. We keep 
the optical power constant, while we take measurements 
for several different detunings A. The radiation pres- 
sure force causes both an optical spring effect resulting 
in a frequency shift of the mechanical resonance, as well 
as damping of the mechanical motion, associated with a 
broadening of its linewidth (see equations (A[4| and (A[5|). 
By fitting the data shown in figures |5^ and (5)3, we can 
extract 71 = 172 Hz and go = 27t x 750 kHz. Knowing 
the mechanical properties of our resonator and the pre- 
cise intracavity photon number, we can now also extract 
the mechanical bath occupancy nb as a function of de- 
tuning from the measured PSDs of the mechanical mode 
(figure [sj); see also [^). This shows us that our mechan- 
ical mode thermalizes to about 16 K for low optical input 
powers, which is close to the cold finger temperature of 
our cryostat of 10 K. 

IV. NOISE SPECTROSCOPY DETAILS 

A. Homodyne measurement with laser noise 

Our experiment is designed to measure the spectral 
density of the fluctuations of the optical field exiting the 
cavity. However any real laser system will have technical 
noise, in addition to the quantum noise associated with 
an ideal coherent source, which adds to the detected noise 
level. Both the signal and local oscillator arm of our setup 
contain this noise which must be taken into account. The 
noise on the signal arm can also be modified non-trivially 
by propagation through the optomechanical system. We 



start by reproducing known results on the operation of an 
ideal, balanced homodyne detection system with signal 
and local oscillator input fields ag and aLo respectively, 
under the influence of noise [28"^ |44l [45] . Most generally, 
these fields consist of coherent tones and OfLO, techni- 
cal (or classical) noise components a^^^it) and aLo,N(^), 
and quantum fluctuations as,vac(^) and aLO,vac(^)- 

as = Qfs + as^N{t) + as,vac(^), (A30) 

ttLO = Q^LO + ^iLO,N(^) + ttLO,vac(^). (A31) 

Since both the local oscillator field and the signal field 
are generated by the same laser, the technical noise on 
the signal and local oscillator will be correlated, and these 
correlations must be accounted for in the analysis. For the 
simplest case, where the signal arm does not experience 
the complex dispersion from interaction with an optome- 
chanical system (e.g. being reflected off the end-mirror far 
detuned from the optical resonator), we expect 

as,N(^) = c^s£,{t) and aLo,N(^) = Q^lo^(^). (A32) 

The function ^{t) is related to the intensity and phase 
fluctuations of the laser light {n{t) and (j){t) respectively): 

a{t) = ao{l ^ n{t))e'^^'^ ^ ao{l ^ n{t) ^ i^{t)) 

^{t) = n{t) + i(j){t) (A33) 

The difference of the photocurrent in the homodyne 
detector is given by 

i{t) = ttstt^Q + aJaLo, (A34) 
which, considering only the technical noise, reduces to 
i{t) = |aLo|#'™'^ + /dc(1 + 2Re{^(t)}), (A35) 

under the assumption that OfLO ^ <^s, using the def- 
initions in equation ( [32| ), and taking the DC current 
/dc = 2Re{a*aLo} = ^jasOfLol cos(^), where is the rel- 
ative phase between the signal and local oscillator. From 
this equation we see that the phase noise (j){t) cannot be 
detected on a balanced homodyne setup. This can be 
understood as being from the detectors fundamental in- 
sensitivity to phase noise on the laser, as the only phase 
reference in the system is the local oscillator, which con- 
tains the same phase fluctuations as the signal. Secondly, 
for the local oscillator phase which makes /dc = inten- 
sity noise is not detected. In a real homodyne detector 
this is only true for a perfect common mode rejection 
ratio (CMRR), which is the case in our setup as the in- 
tensity noise is negligible and the CMRR is > 25 dB. 
For these reasons we use a different setup for character- 
izing the laser phase and intensity noise as described in 
Section ITVBl 

B. Measurement and characterization of laser noise 

In this section we discuss the procedure used for charac- 
terization of our laser (New Focus TLB-6728-P-D). This 
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Figure A5: Optomechanical characterization. We characterize the behavior of the optomechanical system in order to 
extract several parameters such as the intrinsic mechanical linewidth 7i, the optomechanical coupling rate go, and the bath 
temperature Tb (rib), a, The effective mechanical frequency ujm = cjmo -\- Su described in equation (A|4]) is plotted as a function 
of the laser detuning A = cjo — cjl (shown here in units of wavelength) . The frequency shift is due to the optical spring effect 
caused by radiation pressure, b, The optomechanical interaction also causes the intrinsic linewidth 7i of the mechanical mode 
to be broadened as the detuning is changed (cf. equation (^^). c, The area under the mechanical Lorentzian is also modified 
depending on A, and is shown here, normalized to shot-noise. The fits (green lines) in a-c are now used to obtain 71, go and 
rib (see text for details). The plot in d shows a normalized cavity scan, which is used to determine the exact detunings in a-c, 
with every red data point corresponding to a data point in a, b and c. 
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Figure A6: Experimental setup for characterization of 
intensity and phase noise. The laser is amplitude stabilized 
and an attenuator is used to select the desired optical power. 
For the phase noise measurement the light is sent through a 
Mach-Zehnder interferometer (MZI) with a free spectral range 
of 115 MHz. The laser is locked to the center of the interfer- 
ence fringe allowing frequency noise to be converted to inten- 
sity noise. The light is then detected on a New Focus Model 
1811 photodetector and the photocurrent detected on a spec- 
trum analyzer. The same setup is used to detect intensity 
noise without the MZI. 



characterization was done using an independent setup, 
shown in figure [6j and involved two measurements di- 
rectly detecting the light. 

The first measurement is to characterize the intensity 
noise where the laser light is sent directly onto a photode- 
tector with the incident power varied. From the theory 
we expect for the detector photocurrent 



m = 



with 



^(LO,vac) 
=0 

|2 



t- ttLO,vac 

/dc(1 + 



(t))t(h.c.) 
2Re{^(t)}), 



The spectral density of the current is then given by 

Sii{u;) = laLol' (1 + \aLo\^Snn{^)) , (A36) 



where 6'nn(^) is the PSD of the intensity noise fluctua- 
tions. For a real detector, this equation is modified by 
the presence of a dark current 5'dark('^) and non- unity 
efficiency (r^det < 1): 

Sn{uj) = 6'dark(^) 

+ |ceLOp (1 + Vdetl^Lol^ Snn{uj)) • (A37) 

We subtract the dark current (measured with the laser 
turned off) from the data, and set bounds on the mag- 
nitude of the intensity noise present in the laser by ex- 
amining the linear and quadratic dependence of the noise 
floor with respect to power. The linear component is due 
to shot-noise, while the quadratic variance is due to the 



intensity noise fluctuations (see equation (A 37)). The re- 
sults are shown in figure [7^ and c. The noise floor was 
seen to only increase linearly with laser power, confirm- 
ing the absence of intensity noise at the frequencies of 
interest. 

A second measurement is done to characterize the 
phase noise properties of the system. By sending the 
laser through a Mach-Zehnder interferometer (MZI) with 
transfer function I{t) = /o(l + sin(27ra;/cjFSR)), the in- 
tensity of the transmitted light will contain fluctuations 
related to the frequency fluctuations of the light (see fig- 
ure [t]). The free spectral range (FSR) of the MZI is 
^fsr/^tt = 115 MHz. For a real detector, and assum- 
ing uj <C cjfsr, we arrive at 



Let— 2 



Sll{uj) = 5'dark(^) + |Q^Lo|^ ^1 + Vdi 

Some phase noise was detected, as shown in figure [7)3 
and d and the quadratic dependence of the PSD on sig- 
nal power. The spectral densities show a roll-off due to 
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Figure A7: Laser noise characterization, a, We measure the power spectral density (PSD) of our laser for several powers, 
normalize to and then subtract it from the dark current of the detector. The same measurement is performed using a Mach- 
Zehnder interferometer locked at half of the fringe amplitude in order to convert any frequency noise to intensity noise to allow 
detection and is shown in b. c, Plot of the mean value of the PSD around the mechanical frequency cjm from the measurement 
done in a as a function of power. The good linear fit (red line) indicates that no intensity noise is present, d, Mean PSD of 
the measurement in b. The quadratic fit (red line) shows that phase noise is indeed present (see text for more details). 



the FSR of the MZI. It was found that in the frequency 
range of interest, 1 MHz < a;/27r < 40 MHz, the fre- 
quency noise spectral density, Sujuj{^) = uj'^ S (/)(/) (uj)^ is flat, 
and roughly equal to 3 — 6 x 10^ rad^ - Hz, in agreement 
with previous characterization of the same lasers at higher 
frequencies [41^ . 



C. Linearity of detector with local oscillator power 

Having characterized the laser with an independent 
setup, we try to understand the properties of the mea- 
surement system. Our first measurement is designed to 
characterize the linearity of the detector and amplifier. 
With /dc = 0, and no signal in the signal arm of the 
BHD, we expect the system to faithfully reproduce the 
relation ( [35| ) showing a linear relationship between local 
oscillator power and the detected signal vacuum fluctua- 
tion (shot-noise) noise level. It is observed that the mean 
value of the PSDs linearly depend on the input power as 
expected and shown in figure [8) This indicates that our 
detector (and its amplifier) are in fact linear. The red 
line is a linear fit, with a coefficient of determination of 
= 0.9999. Although we already confirmed that no 
measurable amount of intensity noise is present (cf. fig- 
ure [t]) , in the case we would have an appreciable amount 
of noise this measurement would show that it is smaller 
than the CMRR. 



a b 




Frequency (MHz) Local Oscillator Power (mW) 

Figure A8: Noise level versus power, a, Electronic noise 
power spectral densities from the balanced homodyne detector 
at different local oscillator powers (under a balanced condi- 
tion). The red trace corresponds to the electronic noise floor 
with zero local oscillator power, i.e. the dark current, b, Mean 
value of the power spectral densities shown in a as a function 
of local oscillator power. In this plot the electronic noise or 
dark current contribution (0.12 pW/Hz, shown by the dashed 
black line) is subtracted. The red line is a linear fit, which has 
a coefficient of determination = 0.9999. The local oscilla- 
tor power used in the experiment presented in the main text 
corresponds to 3.0 mW. 
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D. Detected noise level with unbalancing 



in a modification on the order of one percent. 
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Figure A9: Amplifier gain, a, Shown are the power spec- 
tral densities (PSDs) of the local oscillator as a function of 
the balancing of the optical power in the two paths of the ho- 
modyne detector. Each trace represents a different ratio of 
power in each path. These traces were taken with a local os- 
cillator power of 3.0 mW, as used in the experiment, b, The 
mean value of the PSDs normalized to the perfectly balanced 
PSD are shown as a function of the difference voltage on the 
two photodiodes in the homodyne detector, where zero volt- 
age represents perfect balancing. The green line is a linear fit 
to the data, while the black curve is a quadratic model, which 
describes any classical intensity noise that could cause the dif- 
ference in the level of the PSDs. The red curve is the sum 
of the two. The change in PSD with homodyne unbalancing 
can be fully explained by the small signal gain weakly depen- 
dent on the detector unbalancing (linear fit) and no classical 
intensity noise (as previously determined). 



A second measurement with vacuum input on the signal 
is done to understand how the amplifier in the homodyne 
detector depends on the DC level of the electronic sig- 
nal after the photocurrent is subtracted. Here we use the 
variable coupler to change the splitting ratio and cause an 
imbalance between the optical power levels in the arms. 
The detected noise floors are shown in figure [9^, and the 
mean detected PSDs are shown in figure |9]3, normalized 
to the shot-noise level. We find that at larger Vdc (I'e- 
lated to /dc linearly), there is a very small (< 2%) drop 
in the gain of the detector. Using a linear fit, we ex- 
tract an adjustment to the gain vs. output DC current. 
This means that for a measured noise power spectral den- 
sity 5'meas(^) taken at a DC voltage Vhc^ we estimate 
that the actual PSD, compensating for modified gain, is 
S{uj) = (1 + FDc/(-0.0096))-^6'meas(c^). This modifica- 
tion is used from here on, and only reduces the amount of 
squeezing we observe, as the quadratures with squeezing 
are always at positive voltages. Additionally, the largest 
DC voltages we work at are roughly ±1 V, which results 
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Figure A 10: Detuned noise. The laser is detuned with re- 
spect to the cavity resonance by 1 nm and spectra are taken 
using the homodyne detector over a range of phase angles, 
with a local oscillator power of 3.0 mW. This lets us estimate 
the amount of additional intensity noise we might acquire in 
our optical signal train, a. The dotted blue line shows the 
amplitude of the interference of the signal and local oscilla- 
tor as a function of time. We lock at several relative phases 
(color-coded from green to red in a and b) and plot the asso- 
ciated normalized power spectral densities (PSD) relative to 
shot-noise in b. For every second measurement we switch the 
signal beam off to obtain the shot-noise level (blue traces). 
The maximum difference in the noise level is around 0.5%. 



E. Estimating added noise in the optical train 

In our third measurement, we reflect the laser light off 
the end mirror of the waveguide coupler (detuned by 1 nm 
from the cavity), and measure the detected noise level as 
a function of ^, the phase difference between local oscil- 
lator and signal. This measurement is sensitive to both 
the conversion of phase noise to intensity noise through 
dispersion in the optical train, and added noise due to ad- 
ditional noise processes in the optical train such as guided 
acoustic- wave Brillouin scattering (GAWBS) [46 , which 
could cause uncorrelated noise in the local oscillator and 



signal arms (see Eq. (A35)). The results of this measure- 
ment are shown in figure 10 i and b. The first figure shows 
the DC interference signal between the local oscillator and 
signal used in the measurement. The LO power used for 
this experiment was the same as for the actual squeezing 
data, and the signal level used is on the same order as 
used for the highest power measurements, as is evident 
from the swing of about 1.4 V in the DC interference sig- 
nal. The highest DC swing observed in the experiment 
was 1.6 V. The second figure shows the normalized (to 
shot-noise) power spectral density where an increase of at 
most 0.5% is observed, indicating these sources of noise 
do not contribute in our experiment. 
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F. The effect of laser phase noise 

Using the measured value for the spectral density of 
phase and frequency fluctuations from section [TV B[ the 
effect of laser technical noise on the detected squeezing 
spectra can be calculated. Following the derivation in 
section I C and taking the classical noise component of the 
fleld input to the cavity to be cl[^\(^) = iaijj(j){uj) (with 

a corresponding LO phase noise of a^Q{(ju) = zaLo^(^)), 
we arrive at an expression for the output noise due to 
input phase noise from the cavity: 



min(l + Ai{uj) - A2{uj)muj). (A38) 



Without optomechanical interaction {G = 0) we flnd 
Ai{uj) = -i^/{i{A - cj) + and A2{uj) = 0. We cal- 

culate the expression for the current noise due to laser 
phase noise using this expression: 



^LO^out (^) 

F{uj)(l){uj) 



«L0[4ut(-^)]* 



(N)/ X , r (N)/ 

out^LO + ^^outi^^LO r 



-UJ 



(A39) 



with F{uj) = i\ai^oain\[e~'^{r{uj) - r(0)) + e'^(r(-cj) - 
r(0))*]. The PSD of the photocurrent due to phase noise 
is found to be 



sf^\uj) = \F{oj)\'S^^{oj). 



(A40) 



For a system with no dispersion, r{uj) = const., it can be 
easily shown that F{uj) = as expected. For an over- 
coupled cavity with no optomechanical coupling, r{uj) = 
1 — hv/{i{A — uj) -\- so r(co') — r(0) ~ Aiuj/n^ and we 

have F{uj) = 8i\ai,oai^\sm{6){uj/ hz). The uj dependence 
of F{(j) means that a flat frequency fluctuation spectrum 
cx cj"^, as we observe) adds a flat noise floor to the 
detected signal. 

Finally we note that phase noise on the laser can drive 
the mechanical motion and cause heating. This effect is 
negligible since we are tuned near resonance, where only 
the intensity fluctuations affect the mechanics, and our 
cavity has a very large linewidth tz. 



G. Phenomenological dispersive noise model: the 
effect of structural damping 



Q = ^(cJmo)"^- We are interested in the wideband prop- 
erties of ^(cj), since the spectral properties of the thermal 
fluctuations are related to the spectrum ^(cj), following 
the fluctuation-dissipation theorem. 

In the case of our experiments, we observed noise floors 
for Sxx following a cj~^ power law on the low frequency 
end. This sort of noise power law corresponds to a flat 
spectrum for the lag angle ^(cj) = const, over the fre- 
quency range of interest. Unlike viscous damping which 
can be simply shown to have ^(cj) oc uj (since the force 
is proportional to velocity), a lag angle constant in fre- 
quency lacks a simple physical explanation, though it is 
ubiquitous in many types of mechanical resonators and 
commonly called "structural damping" [34] . 

In the input-output formalism outlined in section |T] we 
model this type of noise by taking the mechanical damp- 
ing rate 71 to be spectrally flat, and using frequency 
dependent bath correlation functions {bin{^)bl^{^')) = 
{n{uj) + l)d{uj + uj'), {bljuj)kr,{uj')) = fi{uj)S{uj + uj'), 
{bl{uj)bl{uj')) = 0, and {bMknio^')) = 0. This con- 
stitutes our single-mode thermal noise model. 

In any real optomechanical system, a family of mechan- 
ical modes couples to the optical resonance. In the modal 
picture which we use here, each of these mechanical reso- 
nances can be thought to add to the detected noise floor 
with its contribution scaling at the low-frequency end as 
uj~^. The contribution of each mode is proportional to 
the bath temperature, g^^, 7i^k7 and uj^\. We lump 
all of these contributions into a single effective mechani- 
cal resonance, with its properties (not all independent) 
determined by fltting to the low frequency end of the 
noise floor. This mechanical resonator is modeled with 
a mechanical frequency cjm/27r = 50 MHz (so we oper- 
ate in the low frequency tail) , a mechanical quality factor 
Qm = 100, and a total coupling rate of ^o/27r = 100 kHz. 
We found that this model reproduced the magnitude and 
phase (the quadrature in which the noise is detected) of 
the uj~^ noise well, if an additional intracavity photon- 
dependent heating of cq = 3.2 x 10~^ K/photons is as- 
sumed. These background noise floors are plotted in flg- 
This cavity heating rate leads to the effective bath 
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ure 

temperature to nearly double at the highest input powers, 
going from 16 K to over 30 K. This amount of heating is 
in line with what we expect from thin-fllm photonic crys- 
tals we have fabricated in the past operating in the same 
cryostat |24] . 



Mechanical damping of resonators and the associated 
fluctuations from coupling to the thermal bath has long 
been considered as an impediment to measuring weak 
forces in gravitational wave detectors [33l |34l [47H49] . In 
these studies the effect of the bath has often been en- 
capsulated in a parameter ^(cj), representing the lag 
angle in the response of the material to a force. This 
lag angle is the complex part of the spring constant: 
F = —k{l^i'^{uj))x. The quality factor of the res- 
onator is given by the narrow-band properties of the lag 
angle and its value at the mechanical resonance frequency. 



H. Phenomonological absorptive noise model 

In addition to the noise in the quadrature of the me- 
chanical motion (which arises from fluctuations in the 
cavity frequency cJq, and we suspect is mechanical in ori- 
gin), we observed a signiflcant amount of noise in the 
opposite quadrature, which can be interpreted to arise 
from fluctuations of the cavity decay rate hz. Additionally, 
we observed a different noise floor power law (cj"^/^) for 
this noise, which may rule out an optomechanical origin. 
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Figure All: Noise model and experimental results. The 

complete noise model, and constituent components, are plot- 
ted and compared to the experimental, shot-noise-subtracted 
power spectral density (PSD) for a quadrature sensitive to the 
mechanical motion (solid curves) and an insensitive quadra- 
ture (dashed curves). The black lines are the experimental 
PSDs. The red lines represent the full noise model includ- 
ing contributions from a single mechanical mode (blue line), 
phase noise of the laser (brown line), a nd the extra thermal 
noise (green line) as described in section IV G The deviation 
between the modeled and experimental data predominantly 
results from additional mechanical modes. 



Figure A12: Power spectral densities (PSD) of noise 
contributions with varying powers. The complete noise 
model along with its constituent components and experimen- 
tal data are shown for varying optical powers. All curves have 
been normalized to the shot-noise level. The experimental 
data are shown in black with the full noise model in red con- 
sisting of the single mechanical mode (dashed blue), phase 
noise (dashed brown), and extra thermal noise (dashed green). 
The optical power scaling is represented by the transparency 
of the individual curves with curves becoming less transparent 
with increasing optical power. 



The power law scaling agreed with thermorefractive noise 
studied extensively in the context of gravitational wave 
detection [48], microspheres [50 , and microtoroids [51], 
but it is expected that thermorefractive coupling is pre- 
dominantly in the same quadrature as the mechanical 
noise, which is not observed here. Also, if thermorefrac- 
tive, the noise should show strong variation with temper- 
ature through both a quadratic temperature scaling (T^) 
and an extremely steep variation of dn / dT in the temper- 
ature range of 16 K to 30 K which was not observed. 
At this point, we have no noise model to explain the ob- 
served fluctuations, and the origin of this noise will be the 
subject of further investigation to be presented at a later 
time. A phenomenological noise model was instead used, 
where fluctuations in the cavity linewidth proportional 
to the intracavity power with a cj~^/^ noise spectrum are 
assumed. 



I. Comparing measured spectra to theoretically 
predicted spectra 

Our spectrum analyzer (Tektronix RSA3408B) oper- 
ates by taking Fourier transforms of a time domain signal. 
By windowing a short time sample, and calculating its 



energy spectrum, a power spectral density is constructed. 
The size of the window in the time-domain affects the 
resolution bandwidth, and is well known in signal pro- 
cessing, multiplication by a Gaussian window of length r 
is equivalent to convolution of the frequency domain sig- 
nal by a Gaussian with width proportional to r~^. All of 



the measured data, except that presented in section III C 
was taken with a 40 MHz window and 300 kHz resolution 
bandwidth. The spectra contain 501 points spaced by 
80 kHz in the frequency domain. Additionally, for a few 
data sets we took narrow band spectra (down to 100 Hz 
resolution bandwidth) and found that the results agreed 
over the regions where squeezing was observed. The the- 
ory was calculated at 100 times finer resolution than the 
sampled data (with 50,000 points), and was then down- 
sampled after a Gaussian convolution step simulating the 
operation of the spectrum analyzer. This only effects the 
size of the mechanical peak, and has no effect on the 
frequency ranges where we see sub-shot-noise fluctuation 
spectra. For the thermometry data in section HI C[ since 
we are interested in the mechanical line widths and areas, 
the span was always chosen to be the minimum allowable 
by the RSA, which is twice as large as the linewidth. 
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Figure A13: Detected noise power at a given frequency vs. the lock angle. In these plots, a series of traces is shown 
of the detected noise level at a given frequency (with resolution bandwidth of 300 kHz) as a function of the locked phase ^lock- 
The grey points are the measured data points. The solid lines are the results of models detailed in this text, and the dashed 
lines represent different components of noise present in each model. The red line shows the full noise model, containing the 
transduced thermal brownian motion from the studied mode, the noise due to structural damping present in the system, the 
phase noise, and the phenomological out-of-quadrature noise. The green line is for a model considering all the same noise 
contributions, except the phenomological component. A model considering a system without any thermal noise is shown in 
blue. With no thermal force on the mechanical systems, the detected signal in this case can be attributed to radiation pressure 
shot-noise heating. The shot-noise level is denoted by a grey line. The contribution due to thermal motion of the mode of 
interest is shown by the dashed blue line. The noise contributions due to phase noise and structural damping are much smaller 
and shown by the green and red dashed lines, respectively. 



SUMMARY OF NOISE MODEL 



spectra shown in the main text. 



In Table [T| we present a summary of the parameters 
used in the theoretical model for the wideband squeezing 
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Figure A14: Close-up of detected noise power at a given frequency vs. the lock angle. This close-up shows regions 
of squeezing, and the colors are the same as in Figure [13] 



VI. MATHEMATICAL DEFINITIONS 



respectively. The Hermit ian conjugate of operator A{t) 
is given by A'^{t) which has the Fourier transform A'^{uj). 
We present here the notational conventions used This is related to A{uj) 



throughout this work for reference. The Fourier and in- 
verse Fourier transforms of operator A{t) are defined as 



A{t) 



27r 
1 



doj e-*"*i(w) and 

) 

dt e^^'i(t), 



(iM)^ = it(_^). 



(A42) 



In the derivations presented here, we typically express a 
given operator in terms of the "input" bath operators. 
f \A-t\ Expected values are then defined as {A{t)) = Tr[pinA(t)], 
where pin is the mixed state describing the bath. Spectral 
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TABLE I: Model Parameters 



Symbol 


Name 


Value 


Measurement 


n X 1 

^joptical 






-L/ow-powei opuicai spec- 
troscopy with waveme- 
ter. |HI B 




Cavity- waveguide coupling efficiency 


0.55 


Low-power optical spec- 

UiUbUUpy WlUii WcLVclllc Lei . 

Verified phase response with 
VNA to distinguish from 
under-coupling. HIB| 


7i/27r 


Mechanical linewidth 


172 Hz 


Linewidth measurement vs. 

icLocI Lie U U-lillig, 111 UlicI iiiUiii- 

pfrv mppi cinrpmprif i <i;pp t;pp- 

tion HIC). 


Po/27r 


Optomechanical coupling rate 


750 kHz 


Linewidth and mechanical 
frequency measurement vs. 
laser detuning in thermom- 
etry measurement (see sec- 
tion HIC). 


rriO 


Bath temperature 


16 K 


Calibrated areas in ther- 



mometry measurement (see 



section HI C ) 



Co 



Heating by optical absorption 



3.2 X 10"^ K/photon 



Rise of cj"^ noise floor 
with optical power (see sec- 
tion IV G). The cavity tem- 
perature according to this 
model rises from 16 K to 
roughly 30 K at the highest 
powers. 



6 X 10^ rad^Hz 



Frequency noise spectral density 



Frequency noise measure- 
ment with Mach-Zehnder 
Interferometer (see sec- 



tion IV B). 



Laser detuning (red laser is positive) (0.044 zb 0.006) At 



The intensity of the reflected 
light is used to initially set 
the detuning. For a more 
accurate determination, the 
value of A minimizing the 
detected signal for the ob- 
served_^ck is found (see Sec- 
tion IHCI. 



lock angle 



The lock point (as in fig- 
ure |3| is used to find the 
phase angle between the 
light reflected from the cav- 
ity and the local oscillator. 



critical lock angle 



This is the lock angle were 
no mechanical signal is de- 
tected. It is found by looking 
at the area of the mechanical 
mode as a function of ftock 



(see Section II C ) 



densities are found by taking the Fourier transform of the 
auto-correlation 



For a real operator, Saa{^) = Saa{—^)- For non-real 
operators, the spectral density can be symmetrized to 



Saa{^) 



dr e^^^(it(t + r)i(t)). (A43) 
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Saa{uj) = 1{Saa{uj) + Saa{-uj)). (A44) 

These spectral densities can also be calculated from the 
frequency domain operators: 

/oo 
dco' (it(w)i(a;')>. (A45) 
-OO 



We use the convention 

X^J^ =dje-'^ ^d]e'^. (A46) 

to define a measured quadrature of the field. This defini- 
tion corresponds to having a phase difference of between 
input light and local oscillator. 
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